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Abstract 

We investigate the mKdV hierarchy with integral type of source (mKdVHWS), which 
consist of the reduced AKNS eigenvalue problem with r = q and the mKdV hierarchy 
with extra term of the integration of square eigenf unction. First we propose a method to 
find the explicit evolution equation for eigenfunction of the auxiliary linear problems of the 
mKdVHWS. Then we determine the evolution equations of scattering data corresponding to 
the mKdVHWS which allow us to solve the equation in the mKdVHWS by inverse scattering 
transformation. 

Keywords: soliton equation with integral type of source, Lax representation, the inverse 
scattering method. 



1 Introduction 

The soliton equations with integral type of source have important physical application, 
for example, the nonlinear Schrodinger equation with integral type of source (NLSEWS) 
is relevant to some problems of plasma physics and solid state physics It was shown 
in [Q] that the NLSEWS can be integrated by the inverse scattering method for the Dirac 
operator. The key point of the application of the inverse scattering method to integration of 
the NLSEWS in ||] is the use of the determining relations playing the same role as different 
operator representations of the Lax type of nonlinear evolution equations integrable by var- 
ious modifications of this method. Just using the determining relations Mel'nikov obtained 
the evolution equations for all the scattering data of the Dirac operator corresponding to 
NLSEWS. Similar method was used to investigate the KdV equation with integral type of 
source (KdVWS) in ||. The reason for use of the determining relations in [2,3] is that the 
evolution equation of eigenfunction for eigenvalue problem corresponding to the NLSEWS 
and KdVWS was not found. In fact, establish of these determining relations and derivation 
of the evolution equations for all scattering data in [2,3] are quite complicated and required 
some skill. 
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In present paper we investigate the new mKdV hierarchy with integral type of sources 
(mKdVHWS), which consist of the reduced AKNS eigenvalue problem with r = q and the 
mKdV hierarchy with extra term of the integration of square eigenfunction. We first present a 
method to construct the zero-curvature representation for mKdVHWS by finding the explicit 
evolution equation for eigenfunction of the auxiliary linear problrm for mKdVHWS. Then 
we present a way to determine the evolution equation for the scattering data corresponding 
to the mKdVHWS , which implies that the mKdVHWS can be integrated by the inverse 
scattering method. Comparing with the method by using determining relation in [2,3], the 
method proposed in this paper for determining the evolution equation of the scattering data 
is quite natural and simple. This general method can be applied to other (1 + 1) -dimensional 
soliton equations with integral type of source. 

2 The mKdV hierarchy with integral type of source 

Consider the reduced AKNS eigenvalue problem for r = q Q 



M -u( 01 V U=(~ X I ). (2.1) 



2 / „ V V2 



The adjoint representation of (|2.1|) reads || 



Set 



Eq. (O) yields 



and in general 



where 



V x = [U, V} = UV- VU. (2.2) 



i=0 



a = -1, b = c = ai = 0, &i = c x = q 

1 2 h 1 
^2 = -q ,&2 = -02 = ~2 qx ''" ' 



^2m+l — C2m+l — £&2m-lj ^2m — — £2m — —^D^m-l; 



a 2m+ i = 0, a 2m = 2D 1 qb 2m - (2.4) 



L = -D 2 -qD- 1 qD, D = DD' 1 = D^D = 1. 

4 ox 
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Set 

2n+l 

4 A 2n+l-i ^ 



2n+l / x 
y(2n+l) _ [ a i °i ] 



i=0 

and take 

(fcL-^tt)- (2 - 6) 

Then the compatibility conditions of Eqs. ( [2.1|) and (|2.6|) give rise to the mKdV hierarchy 

E 

5H 

%n+i = -2fe 2 n+2 = Db 2n +i = D — -p±-, n = 0,l,---, (2.7) 

where 

2fl2n+2 



H 



2n+l 



2n + 1" 



i 2 -^, L(0 2 -^ 2 )=A 2 (0 2 -^ 2 ). (2i 



Using (|2.1), we have 

The eigenvalue problem (2.1) with q vanishing rapidly as |x| tends to infinity has not discrete 
eigenvalue. As proposed in [2,3,6,7], the mKdV hierarchy with integral type of source is 
defined by 

/OO 
C7(U)($(M,0-$OM,C))dC] (2.9a) 
-oo 

4>i,x = -K4>i + #2, fax = q<t>i + i(4>2 (-00,00) (2.9b) 

we assume q(x, £2n+i) tends rather quickly to zero as x — > ±00. According to this condition 
we assume that 

4>i(x,t,Q ~ a(t, ()exp(— i(x), (j) 2 (x,t, Q ~ b(tX)^xp(i(x), x — > —00 (2-10) 

where C = C(t,(), a=a(t, Q and b—b(t, () are complex functions of £ > and £ G (—00, 00). 
Moreover we assume that the functions C, a and b are chosen so that the right-hand side 
of equation (|2.9| ) determines the function absolutely integrable over x along the whole real 
axis. One can easily verify that the requirement will certainly be satisfied if the function E 
and T of the form as argued in 

E = \C(t,0\[\a(t,0\ + \b(t,0\} 2 

r = \^[C(tx)a\tx)}\ + \-^iC(t,Ob 2 (t,0}\ 

at any £ > satisfy the condition 

/oo 
[£(U) + r(U) + r 2 (U)R<oo 
-00 
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3 The Lax representation 

Following the method proposed in [6,7,8], in order to find the zero-curvature representation 
for (|2.9| ), we first consider 

POO 

D[b 2n+1 + / C(t,O(<t>l(x,t,Q-<%(x,t,O)dC] = (3-la) 



k,x = -*C0i + #2, 4>2,x = q<j>i + i(4>2 (-00,00) (3.1b) 



We can obtain the Lax representation for ( |3.1|) by using the adjoint representation ( 2.2 
According to Q2.4|) , ( |2.8|) and (|3.1| ), we may define 



d{ (l{ 7 bi bi^Ci Cj, I 0, 1, • • • , 27Z, 
poo 

b 2 n+2m+l = C 2 „+2m+l = Lb 2n +2m-l = ~ {i() 2m+2 C{t, £) [4>\ {x,t, () - (j)\(x, t, C)]d( 



00 

00 



&2n+2m+2 = -C 2n+2m+2 = ~]:Db 2n+2m+1 = ~ / «) 2m+1 C(t, C) [0?(x, t, C) + ^OMjOMC 



Then 



2 

a 2 n +2m+2 = 2D- 1 ^ 2n+2m+2 = / «) 2m+1 C(t, ^^(x, t, C)0 2 (x, t, ()d( 

J —00 

0-2n+2m+l =0, m = 0, 1, • • • , 

M (2n + i) _ ( B^) \ _ 2n+1 ^f~a k b k \ h (9 \ 

7V ~ ^ C^ +1 ) _D (2n+1) J =A ^c fc -aJ A + 1,0 9j 

where 9 is some constant and 

/1 (2„+l) ,2„+l-l , fl , r° gliCKjgjgfeWUf. t, gjfefa t, 1?) 

= & (O 1 -(*!)" ' 

R (2n+1) V^a ^2n+l-fc^ /"^ <« ~ ^)C(t, ^jfe g) ~ <« + M^gft, ^(g, t, 7]) 

B -}^b k \ + _ 



fc=0 
2n 



«) 2 - (^) 5 
W2ra+1) _ \ „ \2n+l-k , 



fc=0 

2n 



n (2„+i) V- ,, x2 „+i-fc , fl / 2«)(i7y)C(t, 77)01(2:, t, ?y)0 2 (g, t, 77) 



k=0 



k) 2 - (^y 

«) 2 - (^) 2 
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also satisfies the adjoint representation (|2.2|), i.e. 

Ni 2n+ V = [U,N {2n+ % (3.2) 
which, in fact, gives rise to the Lax representation of ( |3.1| ). Since ( |3.1| ) is the stationary 



equation of (|2.9|) , it is easy to find that the zero-curvature representation for the mKdV 



hierarchy with integral type of source ( |2.9| ) is given by 

Ut 2n+1 ~ + [U, N^] = 0, (3.3) 

with the auxiliary linear problems 



1>i\ ( -A q \ \ = ( -K q \f V>i 

i>2 L \q a U Kg K A fa 



(3.4a) 



where A = iQ and 



k=2n 



(x, Wi , C) = (^ (2 " +1) + 0)^i + B {2n+1 V2 = Mi + ^2) A 2 " +1 ~ fc + 0Vi 



fc=0 



OG 



iCC(t v) 

/ ]jQ2~(~y [ 2 ("7)<M X > »7)02(x, i, ?7)^i + («C— *, ^)V ; 2-«+^)0?(^, v)fa\ d V, 



oo 

fc=2n 

_ ./. _ ./. \\2n+l— fc i i-'^ 

OO 



(x, Wi, = C( 2n+1 Vi + (-^ 2 " +1 ) + 9)^ 2 = Mi ~ a k ^j 2 )\ 2n+l - k + Ofo 



-2(i77)^i(x, £, rfifafa v)^2}dr]- (3.4b) 

In this way we find the explicit evolution equations of eigenfunction if). Indeed, this kind of 
evolution equation of eigenfunction was not obtained in [2,3]. 



4 Evolution equation for the reflection coefficients 

Now we can derive equations describing the evolution in time t of the S-matrix elements. 
This can be made as follows. We define the eigenfunctions f~(x, £) = (fi(x, £), f 2 (x, C)) T > 

t(x,C) = (ft(x,0J 2 -(x,C)) T , f + (x,C) = (fHx,0, f2 + (x,C)) T and f+(x,() = (/^,C), 
f 2 {x, Q) T (here and hereafter the "T" means transposition) for the equation (|3.4a|) , and the 
following asymptotics are fulfilled at any ( G (— oo, oo) 

f O*,0 ~ ( J ) e ~^'' r( x 'C) ~ ( °_i ) e ^> os ^ ^ -oo (4.1a) 
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f + (x, ~ ^ J J e**, /+(x, C)~(Jj e"**, as x - +00 (4.1b) 

As is known, the functions f~(x,Q an d f + {x,C) admit an analytical continuation in the 
parameter ( into the upper half-plane Im£ > 0, and the functions f~(x,Q an d / + (x,C) 
admit an analytical continuation in the parameter ( into the lower half-plane lm( < 0. It is 
easily seen that at any real £ G (—00,00) the pair of functions f~(x, £) and f~{x,() forms 
a fundamental system of solutions to (|3.4a| ). Hence, we may define 

f + (x, C) = S 12 (C)f-(x, + S 22 (0f-(x, C) (4.2a) 

; + (x,o = 5 n (or(x,c) + s 21 (or(x,() (4.2b) 

where the quantities S u = S n (C), S 12 = Si 2 ((), S 21 = S 2 i(C) an d 5*22 = S 22 (() are in- 
dependent of x. Taking account of Q4.1Q and (|4.2j ) we get that at any £ G (—00,00) the 
equality 

Sn (0^22(0 -s 12 (Qs 21 (C) = l. (4.3) 

Under the assumption that q(x,t) vanishes rapidly as |x| — > 00, we have 

a = -1, b = c = 0, lim^^aj = lim^^bj = lim\ x \_^Cj = 0, j = 1, 2, • ■ ■ , 2n. 

We denote the parameter # in ( |3.4b| ) corresponding to f + (x,Q by # + and f + (x,Q by # + , 
respectively. Substituting / + (x,C), f + (%,() ( |3.4b| ) respectively, we have 

OT 2n+l fc=0 J-00 V ~ S ?7 + S 

-7r«)H(C) + 7r(iC)^(-C)}A + (^C) 

/*00 J- poo > 

+{J2^0 2U+1 ' k + f -zjH^dn+f ^-H 2 (rj)dr, 

f,_Q J -co V S> J -co V ' S 

+7r(iC)^i(C) -t(O#2(-C)}/ 2 + 0M.C), (4.4a) 



( ^ ( ' ' ^ 5 ^ {f>«) 2 " +1 -* - r -$-rH 2 { V )d V - f -^-MMdri 



0*2n+l , : „ 



2/i 



-7r(iC)^ 2 (C) + 7r«)^i(-C)}/i + (x, C) 



fc=0 



/ —* 1 , iQ — vq iQ + 17] 
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+7T«)# (C) " 7r«)# (-C)}/ 2 + (*, 0, (4.4b) 
+n(zC)H(0 - 7r(t0H(-()}f+(x, 



/CO J- poo J- 

^ Hl (rj)dr,+ f ^H 2 (rj)d V 
-oo V S J —oo 1 + 1, 

-7T(iC)^l(C) + TT^O^C-C)}^^, 0, (4.4c) 

^ = * - £ - £ ^ w 

+7r«)# 2 (C) " TriiQH.i-Oyf+ix, o 

+{ | ( - aiK )---£ ( ^-^_ )WWd , +s - + 



-7r«)//(C) + n(zOH(-()}f+(x, 0, (4.4d) 

where the integral § is taken as the principal value, the quantities 9 + , 6* + will be determined 
in the next section, and 

H (v) = C(t,r])cl) 1 (x,t,r])(p 2 (x,t,r]), 

H 1 (r ] ) = C(t, V ) ( j ) i(x,t, V ), H 2 ( V ) = C(t, V )(Pl(x,t, V ). (4.5) 
As x — > — oo, we find that the following asymptotics are valid: 

f -^jH^dv ~ 7r(^C)C(t,C)a 2 (C,t)e- 2 ^, 

J —oo V S 

/°° -^-MMdri ~ 7r(<C)C(t, -C)a 2 (-C, t)e 2 ^, 

/°° -^-H^drj ~ -7r(i0C(f , C)& 2 (C, t)e 2 ^, (4.6) 
/°° -^-H^dr) ~ -^C)C(t,-C)& 2 (-C,t)e- 

/-oo 1 + C 



Substituting ( [4.2| ) into ( |4.4| ) and using (f4.6|), as x — > — oo, we have 



dt 2n +i J-oo V~C V + C 

-n(z()h(0 + n(zOh(-0}S 22 (0 
-{2 7 r«)/i 1 (C) - 2n(i()h 2 (-0}S 12 (0, 

§^ = {(iC) 2n+1 + e++f(^- + ^-)h( v )d v 
dt 2n+1 77 — C v + C 

+7r (i()h(() -7r(i()h(-()}S 12 (0, (4.7) 

§^ = {-«) 2n+i + ^ + - r (-^ + 4?)^ 

+7r (iC)/i(C)-7r(iC)M-C)}5 2 i(C) ; 
^11 (0 _ {(<)2 n + i + ^+ + /"(_<_ + -^-)h( V )dr, 



where 



dt 2n+ i /-oc ^-C »7 + C 

-ttCiOMO + ^OM-C)}^^) 

-{27r(tC)^(C) " 2 7 r«)/i 1 (-C)}^i(C), 

h(rj) = C(t, V )a(rj,t)b(ri,t) 
h 1 (r,) = C(t,r ] )a 2 ( V ,t), h 2 (r ] ) = C(t, V )b 2 (r ] ,t) 



One can easily see that ifC = 0ora = 6 = then the resultant system (|4.7|) coincides with 
those equations which appear in the case of the mKdV hierarchy without a source. Also one 
can verify that system ( f4.7|) is consistent with equality (f4.3|) . Using (|4.7| ), we find that the 
reflection coefficients 



satisfies the equation 



Ot 2n+ i 7} -Q V + C 

-7r(<0/»(C)+T(iC)^(-C)}i2i(0 

-{27r(zC)/i 2 (0 - 27r(iC)/ii(-C)}, (4.9a) 
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= 2{-(z C ) 2 " +1 - fi-^ + -^Mrfidr, 

ot 2n +i J-oo v~C v + C 

-7r(<C)/»(C)+T(tOfc(-C)}i2i«(0 

-{27r(iC)/ii(C) - 27r«)/ i2 (-C)}. (4.9b) 

Then, it follows from ( f4.9|) that the evolution of the reflection coefficients R±,R2 are in- 
fluenced by the integral type of source which is integration of the square eigenf unctions 
belonging to the continuous spectrum of the spectral problem ([2.1|) . For the case r = q, 
there is no discrete eigenvalue for the spectral problem ( |2.1|) if the potential q = q(x, t) tends 
rather quickly to zero as |x| — > oo. The evolution equations for the reflection coefficients are 
presented by (|4.9| ) which implies that the mKdV hierarchy with integral type of source can 
be solved by the inverse scattering method. 



5 Conclusion 

By means of the reduced AKNS eigenvalue problem with r = q which has no discrete 
eigenvalue, we construct the mKdV hierarchy with integral type of source (mKdVHWS). We 
propose a method to find the evolution equation of eigenfunction corresponding to mKd- 
VHWS and further to determine the evolution equation for scattering data which enable us 
to solve mKdVHWS by Inverse Scattering Transformation. Comparing with the method for 
determining the evolution equation for scattering data in [2,3], our approach is quite natural 
and simple. 

It should be noted that the reduced AKNS spectral problem for r = — q may have the 
discrete eigenvalue. In this case, the right-hand side of equation Q2.9a| ) need to be added by 



the sum of square eigenf unctions of (|2.9b[) corresponding to the discrete eigenvalue. We will 
show in the further coming paper that the mKdV hierarchy with these two kinds of sources 
can also be integrated by Inverse Scattering Transformation. 
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